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t*^ ■ Abstract 

. ■ We study the time-dependent scattering of a quantum mechanical wave packet 

\l ' at a barrier for energies larger than the barrier height, in the semi-classical regime. 

>--^ , More precisely, we are interested in the leading order of the exponentially small 

scattered part of the wave packet in the semiclassical parameter when the energy 
density of the incident wave is sharply peaked around some value. We prove that 
this reflected part has, to leading order, a Gaussian shape centered on the classical 
trajectory for all times soon after its birth time. We give explicit formulas and 
rigorous error bounds for the reflected wave for all of these times. 
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1 Introduction 

We consider the problem of potential scattering for a quantum particle in one dimen- 
sion in the semiclassical limit. Let U : M ^ E be a bounded analytic potential function 
such that lim|j,|^o2 ^i^) = 0- Then a classical particle approaching the potential with 
total energy smaller than max^jgR V{x) is completely reflected. It is well known that for 
a quantum particle tunneling can occur, i.e. parts of the wave function can penetrate 
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the potential barrier. Similarly, for energies larger than the barrier height maxx^B.V{x) 
no reflection occurs for classical particles, but parts of the wave function of a quantum 
particle might be reflected. In the semiclassical limit the non-classical tunneling probabil- 
ities respectively the non-classical reflection probabilities are exponentially small in the 
semiclassical parameter. While the computation of the semiclassical tunneling probability 
can be found in standard textbooks like |LaLij . the problem of determining the reflec- 
tion coefficient in the case of above barrier scattering is much more difficult [LaLii IFrFrj . 
In this paper we are not only interested in the scattering limit but also in the questions 
where, when and how the reflected piece of the wave function emerges in an above barrier 
scattering situation. 

The time-dependent Schrodinger equation for a quantum particle moving in one 
dimension in the potential V is 

where < e ^ 1 is the semiclassical parameter. A solution to (H]) can be given in the 
form of a wave packet: for any fixed energy E > max3;gR V{x), let 0(x, E) be a solution 
to the stationary Schrodinger equation 

-e^dl<j>={E-V{x))cf>=:p'{x,E)cl>, (2) 

where 

pix,E) = yjE - V{x) > 

is the classical momentum at energy E. Then the function 

*(x,t,£):= I Q{E,e)e-'*'^^'(j){x,E) dE, (3) 

for some regular enough energy density Q{E,e) is a solution of ^. We will require 
boundary conditions on the solution of ^ that lead to a wave packet that is incoming 
from X — -l-cxD and we will assume the energy density is sharply peaked around a value 
Eq > maxa;gRF(a;). 

In order to understand the emergence of exponentially small reflected waves from ^ 
we basically need to solve two problems. Roughly speaking, we need to flrst determine the 
solutions to ^ with appropriate boundary conditions up to errors sufficiently small as 
e ^ 0. Since we are interested in the solution of ((T|) for all times, we need to determine the 
solution of ([2]) and also need to split off the piece corresponding to a reflected wave for all 
X e R and not only asymptotically for |a;| — > cx). The large |a;|s regime of ([2|) is rigorously 
studied in |J oPf2[ I Raj by means of complex WKB methods and yields the exponentially 
small leading order of the reflected piece. On the other hand, the behaviour for all a;s 
of ([2]) has been investigated at a theoretical physics level in [Bel] , where error function 
behaviour of the reflected wave piece has first been derived. However, our method is quite 
different and allows for a rigourous treatment, which is one main new contribution of the 
present paper. 

In a second step we have to evaluate (|3]) using the approximate (j)s from step one 
and extract the leading order expression for the reflected part of the wave, for a suitable 
energy density Q{E,e). Since we know that the object of interest, namely the reflected 



wave, is exponentially small in e, both steps have to be done with great care and only 
errors smaller than the exponentially small leading order quantity are allowed. This has 
so far been done only in the scattering regime for similar problems |HaJo31 IJoMa] . By 
contrast, we will be able to carry out the analysis for all times just after the birth time of 
the reflected wave, not only in an asymptotic regime of large positions and large times. 
This is the second main new aspect of our work. 

Before further explaining the result, let us give some more details on how to ap- 
proach the two problems mentioned. For solving ^ we use techniques developed in 
[BeTeli IBeTe2] . cf. also |HaJo2j in the context of adiabatic transition histories. The 
connection to adiabatic theory is easily seen by writing the second order ODE @ as a 
system of first order ODEs. Put 






where the prime denotes one derivative with respect to x, then (f> solves ([2|) if and only 
if ijj solves 

i£5.V(a;,i?)- (^^2(f^^) l^^{x,E)=:H{x,E)^{x,E). (4) 

In |BeTe2j we studied the solutions of the adiabatic problem 

iedt^it) = H{t)i;{t) , (5) 

where H(t) is a time-dependent real symmetric 2 x 2-matrix. So the two differences 
between ^ and ([5|) are that H{x, E) is not symmetric, but still has two real eigenvalues, 
and that an additional parameter dependence on E occurs. In Section 2 we explain how 
the results from [BeTe2) on (O translate to (|3]). 

From this we find in Section 3 that the solution to ^ can be written as 

0(a;, E) = 0ioft(x, E) + bright (2:, E) , (6) 

where <j)\cit{x^ E) corresponds to a wave traveling to the left and is supported on all of R, 
while bright (a;, E) is the reflected part traveling to the right and is essentially supported 
to the right of the potential. The difference between ^ and the usual WKB splitting 

e7/o"p('-.-E)d'- e-iIoP'-^^E)dr 

</.(x, E) = /irft(2:, E) ^_^ + /Hght(a;, E) ^==^ (7) 

^/p{x, E) ^/p{x, E) 

is that the reflected part (/)right(a;, E) is exponentially small for all a; € M and not only for 
X — > -l-oo as fiig\it{x,E). Indeed, we will show that 4>^ig\^t{x,E) is of order e~^'=(-^)/^ for 
an energy dependent decay rate £,c{E). In the scattering limit |a;| — > 00 the splittings @ 
and ([7]) agree; so the main point about ^ is that we can speak about the exponentially 
small reflected wave also at finite values of x. 

Inserting ([6]) into ([3]) allows us to also split the time-dependent wave packet 

^(a;,f)= /"Q(^,e)e"^*^/^(0ieft(x,S) + (/.right(2:,S))dS=:*icft(2:,t) + *,ight(2:,t). (8) 



Since bright (a;, S) is essentially of order e-^^^^)/^^ the leading order contribution to 
^right(a;,t) conies from energies E near the maximum E*{e) of |Q(i?,e)|e~^=^^^/'^. In 
order to see nice asymptotics in the semiclassical limit e — > we impose conditions on 
Q{E,e) that guarantee that E*{e) has a limit E* as e ^ 0. Typically, this will be true 
for a density Q{E,e) which is sharply peaked at some value £'0 as e ^ 0, in accordance 
with the traditional picture of a semiclassical wave packet. In particular, it should be 
noted that in general, the critical energy E* is different from Eq, the energy on which 
the density concentrates. 

The main results of our paper are increasingly explicit formulas for the leading order 
exponentially small reflected wave bright (a;, t) not only in the scattering limit, covered by 
[Ha Jo3[ fJoMaj . but for most finite times and positions. In particular these formulas show 
where, when and how the non-classical reflected wave emerges. Let us note that similar 
questions can be asked for more general dispersive evolution equations, in the same 
spirit as the systems considered in [JoMa] . However, the missing piece of information 
that forbids us to deal with such systems is an equivalent of the analysis performed in 
|BeTelirBeTe2j . and adapted here to the scattering setup, that yields the exponentially 
small leading order of the solutions to ^, for all xs. 

Our methods and results belong to the realm of semiclassical analysis, in particular, 
to exponential asymptotics and the specific problem we study has been considered in 
numerous works. Rather than attempting to review the whole litterature relevant to the 
topic, we want to point out the differences with respect to the results directly related to 
the problem at hand. Exponentially small bounds can be obtained quite generally for the 
type of problem we are dealing with, see e.g. [Fell lFe2l [JKH [Si} |JoFf7| iNel iMal iHa^ 
Moreover, existing results on the exponentially small leading order usually are either 
obtained for the scattering regime only, Fcl| IFe2l |JoPf2[ IRa| [Jol IHaJo3|. IJoMaj , or non- 
rigorously [Bell IBe2| IBeLij , or both [WiMoj . To our knowledge, the only exceptions so 
far for the time independent case are [BeTeli [BeTe21 lH"aJo2) . which we follow and adapt, 
but as mentioned above we need to be even more careful here in order to get uniform 
error terms in the energy variable. Finally, we focus on getting explicit formulae for the 
leading term of the reflected wave, rather than proving structural theorems and general 
statements, which justifies our choice of energy density. 

The detailed statements of our results are too involved to sum them up in the 
introduction; let us instead mention a few qualitative features. In the semiclassical limit 
e — > the reflected part of the wave is localized in a -ye-neighborhood of a classical 
trajectory qt starting at a well defined transition time at a well defined position q* 
with velocity ^E* — V{q*). A priori we distinguish three regimes. A y^-neighborhood 
of the transition point q* , which is located where a certain Stokes line in the complex 
plane crosses the real axis, is the birth region, where the reflected part of the wave 
emerges. It remains ^/e-localized near the trajectory qt of a classical particle in the 
potential V with energy E* > maxV^ for all finite times, which defines our second regime. 
Therefore limt_>oo Qt — 00, i.e. the trajectory belongs to a scattering state and moves to 
spatial infinity. The third regime is thus the scattering region, where the wave packet 
moves freely. We give precise characterizations of the reflected wave in all three regimes. 
Theorem [5] It turns out that as soon as the reflected wave leaves the birth region, it is 
well approximated by a Gaussian wave packet centered at the classical trajectory qt, see 
Theorem [HI 



We now give a plan of the paper. In order to precisely formulate our main results, 
we need to discuss in detail under which conditions and in which sense we can solve (|4]) 
up to exponentially small errors. This is discussed in Section 2 where the main result is 
Theorem 131 The proofs are based on the techniques developed in [BeTeli lBeTe21 lB"eTe3j 
and can be found in Sections 5 and 6. In Section 3 we explain how to arrive from our 
solutions of ^ at the main theorems, Theorems I3HS1 by first obtaining solutions of ^ 
split according to ([6]) and then evaluating ([8]). The mathematical details on the asymp- 
totic analysis of (O are covered in Section 4, where the main tool is Laplace's method. 
In Section 5 we explain how to bring ^ into an almost diagonal form, with off-diagonal 
elements of order e" for arbitrary n G N, see Theorem [T] We proceed analogous to the 
well known transformation to super-adiabatic representations in adiabatic theory, see 
[BeTeli IBeTe2] for references. In Section 6 we analyze the asymptotic behavior of the 
off-diagonal terms in the super-adiabatic representations. By optimal truncation we find 
the optimal super-adiabatic representation in which the off-diagonal terms are expo- 
nentially small and explicit at leading order. Then standard perturbation theory yields 
Theorem |31 

2 Superadiabatic representations for non-selfadjoint 
Hamiltonans 

In this section we adapt the results of [BeTel] and [BeTe2) to fit our needs. The main 
difference is that the Hamiltonian matrix we will have to deal with is not self-adjoint; 
this will require some adjustments which we carry out in the Section [6l but fortunately 
much of the theory from [BeTeli IBeTe2] carries over. The equations we consider are 

ied.Me, X, E) = H{x, E)4>{e, x, E) := ( ^^^^^ ^^ J \ ^{e, x, E) (9) 



where p{x, E) = ^E — V{x\ and V is analytic in a neighbourhood of the real axis with 
sup3.gR ^(^) ^ ^- We will here give the optimal superadiabatic representation for this 
system, which is an e-dependent basis transform that diagonalizes ^ up to exponentially 
small errors while still agreeing with the eigenbasis of H{x,E) for |a;| — > cx). In this 
basis, for suitable initial conditions the second component of the solution is given by an 
exponentially small error function with superimposed oscillations, to leading order; this 
fact will enable us to describe the time development of the exponentially small reflected 
wave. 

We start by introducing the invertible transformation 

a^,E)^2 f p{y,E)dy. (10) 

Jo 

Equation ([9]) then transforms into 

ied^i^ie, e, E) = Hit EWe, ^ i?) := ^ ( ^^^ j^^ ^ ) V^e, C, E), (11) 



where p{^{x,E),E) = p{x,E). Let us for the moment consider fixed E and drop the 
dependence on E from the notation. The matrix on the right hand side of pT|) has eigen- 
values — 1/2 and 1/2, and is diagonahzed by the basis transform implemented through 



rtn - ( vW) i/vW) ^ 



We then find 



0^Tiied^-HiO)T-\T^)^(^ied^-l(^ .J^^^ ^"^^ )) ^a, (12) 

where 0'{^) — ^^4^ is the adiabatic coupling function. The following theorem gives the 
n-th superadiabatic representation for every n G N. 

Theorem 1. Assume 9'{£,) G C°°(M). Then for each n e No there exists an invertible 
matrix Tn — Tn{e,£,) such that 

TAied,-H)T-^=ied,-[ ^ (13) 

where P7i{e,^) = ^ + C(e^), kn{e,£,) = {xn{0 ~ ^n('?))(l + ^(s)); '^^'^ ^he Xj, Zj can be 
calculated from the recursive equations 

XI = ^9\ yi = zi = 0, (14) 

(-'n+l = -12„ -t- it7 (/„ (15) 

2j^n+l-j) (16) 
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(17) 



We will give the proof of this theorem in Section [5l There we will even give explicit 
expressions of T„, p„ and fc„ in terms of the x^, yi and z^, but these are too messy to 
write down here. Note that none of them is a polynomial in e. 

As it stands p^ is not useful, since neither the n-dependence nor the ^-dependence 
of the off-diagonal terms is controlled. In order to achieve this, we need to make stronger 
assumptions on 6' . As the present situation is very close to the superadiabatic transition 
histories considered in [BeTe2) . the conditions will be identical: we require that 0' is 
analytic in a strip around the real line, and that the singularities which limit this strip are 
first order poles plus lower order corrections. As explained in [BeTe2) . those corrections 
are controlled using the following norm on functions of a real variable: 

Definition 1. Let ^c > 0, a > and / C M be an interval. For / e C°°{I) we define 

||/||(,^„_^^) ■.= sn^snv\d^f{t)\-^^<^ (18) 

^ '^ ' tei fe>o ' 1 (a-l- k) 

and 



Fa..iM) = {f (^ C^il) ■■ \\f\\il,a..^^) < ^} ■ 



We now state the assumptions that we will make on the adiabatic coupling function. 
The first one will be required for regions of R that are far away from the first order poles 
that are nearest to the real axis, while the second one will be required close to such poles. 
The quantities ^c and ^j. will be the imaginary and real part of those poles. 



Assumption 1: For a compact interval I and k > ^c > let 9' {^) E -Fi.k(/). 
Assumption 2: For 7, fr, ^c G R let 

be the sum of two complex conjugate first order poles located at ^r ± i^c with residues 7. 
On a compact interval I G [^r — Cc: ■Cr + Cc] with ^r € I we assume that 

o'{0 = o'o{0 + ()r{0 with e',{O^Fa,Ui) (19) 

for some 7, fc,Cr € R, < a < 1. 

Under these assumptions, the following theorem holds: 

Theorem 2. Let J^^ C (0, 00), Ja C (0,1) and J-y C (0,cx)) be compact intervals, and 
for given ^c > Ze< < (Te < 2 6e such that 



Tie — 1 + Cg 



is an even integer. 



(20) 



(i) There exists Eq > and a locally bounded function (j)i : R"*" — s- R"*" with (pi (x) = 
0{x) as X —^ 0, such that the following holds: uniformly all in e € (0,eo]) (^^^ 
compact intervals / C R and all matrices H{S,) as given in ilS^) and satisfying 
Assumption 1 for some ^c S J(,^, the elements of the superadiabatic Hamiltonian as 
in lilSp and the transformation Tn^ — T„^ (e, ^) with n^ as in li2U\) satisfy 



Pnj£,0- 2 



< ^'<^i(ll^'ll(/,i„ 



and 



|T„ 



Til < e. 



</>i(ll^'ll(/,iX 



(21) 
(22) 

(23) 



(ii) Define 



3(='^) = 2iJ^sin(f^)e 



e e 2e5c 



(^ 






+ 






(24) 



There exists eg > and a locally bounded function 02 : R^ -^ R^ with 02 (x) = 
0{x) as X —^ 0, such that the following holds: uniformly in all e < Eq, all compact 



intervals J C M, and all H{^) as in il2\) satisfying Assumption 2 for some ^c G J^^, 
a G Ja, J ^ J-y, we have 

|e"^+'fc„.(e,e)-5(e,0| <£'-"e"^</.2(M), (25) 

where M = max-^ ll^'ll(/ if); ll^rll{/ q f ) f • Furthermore, \21]) and \23\) hold with 

In words, Theorem [2] states that we can find a transformation T„ which is e-close 
to the one that diagonalizes H{£^), but which diagonalizes the full equation P^ up to 
exponentially small errors. What is more, for those ^ where the off-diagonal elements are 
close to maximal, they are given by explicitly known Gaussian functions with superim- 
posed oscillations to leading order. The dependence on E, which we dropped earlier from 
our notation, comes in only through 9' , which means through ^c, ^r and 7 in the leading 
term, and through H^'ILj i f \ in the corrections. Since later we will have to consider a 
whole range of values for E, it is of crucial importance that all our error terms are small 
uniformly in basically all parameters, as stated above in Theorem [S] On the other hand 
note that in those regions that are more than ^/e away from ^r, the error term dominates 
in ((25)) . We will give a proof of those parts of Theorem [2] that are not identical to the 
situation in ' BeTeli rBeTe2] in Section [Sj 

Theorem [5] shows how to diagonalize equation (jlip up to exponentially small er- 
rors. We will actually need an exponentially accurate solution to this equation in the 
whole ^ domain. Here we will give a simple version of the corresponding result which is 
analogous to the one in BeTe3J , and for which only mild integrability conditions on 9' 
and its derivatives are necessary. In [BeTe3| it is shown that a sufficient condition is the 
integrability of ^ i-^ ll^'llff 1 £ +5) fo^' some S > 0, outside of a compact interval, and that 
a sufficient condition for this is the integrability of 

^^snp{\9'{z)\:zeC,\^-z\<^, + d} 

outside of a compact interval. 

Theorem 3. Assume that 9' has two simple poles at ±i^c fl^ leading order, and no 
further singularity of distance less or equal to ^c to the real axis. In addition assume 
that for some 5 > Q, S, ^—f ||^'||/£ i^ ,g\ is integrahle outside of a compact interval. Let 
■'/'«(£, — (''/'+,«(£, 0; V'-,n(£, 0) ^s the solution of the equation il2\) in the optimal 
superadiabatic representation, i.e. 

T„^ {ied^ - H) T- V„ (£,0 = 0, (26) 

subject to the boundary conditions 

lim V+,n(£)Oe^ = 1, 1™ V-,«(£,C) = 0. 

Then 

V'-,„(£,0 = 2sin (^) e-^ei-^erf (^^) + 0(£^-"e-^) , 



where 

erf(a;) = —= I e"^ ds 



^"".'-oo 



is the error function normalized to be at — cxd. 

The proof of Theorem [3] requires no changes as compared to the one given in [BeTeSj 
and is therefore omitted. Let us remark already here, however, that for our apphcations 
to the reflected wave we will need precise control over the ^-dependence of the error term 
above, which means we will basically have to reprove a refined version of Theorem [3] 
later. Also for later purposes, we note that Theorem [3] also gives exponentially accurate 
information about ip+y. defining 

wc find using ([26]) (given explicitly in (j38|l ) that 



\i^+AC} - i^+MO\ 



-f e"+ifc„,(e,y)V^„,„(£,2/)dy =o(£-i/2e-^). (27) 



The last equality follows from (P^ and ([M)l and the integrability condition in Theorem[3J 
We end this section by linking the conditions on 6' to conditions on the function p{x, E) 
appearing in ((9]). The first condition is 

(Al) We assume that p{x, E) > c for some c > and all a; G M, and that p(-, E) is 
analytic in a neighbourhood U of the real axis. 



While the map x i-^ p{x,E) is analytic on U by (Al), it may have singular points 
in the complex plane. For each such point zo we define a straight half line Bz^ pointing 
away from the real axis. In the case that a branch cut of p is needed at zq we take it to 
lie on Bzo- Since Uq := <C\ [J^ B^a is simply connected, the function 

az,E)^2 f piy,E)dy, (28) 

Jo 

for z £ Uq is single- valued if we require that the path of integration has to lie entirely in 
Uq. ^{z, E) is the unique analytic extension of ^ as given in (fTU]) to Uq. For a; G M we put 

CAx):^{zeUo:\a^)~ax)\<r}, 

denote by Cr.oix) the connected component of C,(a;) containing x and put 



R{x) := sup{r > : Crfl{x) C Uq}, 
and Rq — infa;gRi?(x). Finally we put 

Uoo = U Cr,^o{x). (29) 



xei 
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The point of this construction is that now the set dUoo \ [/q contains only those points 
on the boundary of Uq that are closest to the real axis in the metric given by 

d(z,z') = ICW-e(^')l- (30) 

It turns out that these are the points on which we need to impose conditions. 



(A2) We assume that dUoo\Uo — {zcrit, ^^crit} with Imzcrit > 0, and that for z G C/qo 
near Zcrit we have 

p(z, E) = iz- Zcrit{E)f^''^ {A{E) + g{z - Zerit, E)) for some PiE) > -1, f3{E) ^ 0, 

where g is analytic near z = with 17(0, E) = and A{E) ^ 0. Furthermore, we assume 
that there exists J > such that for all other critical points zq of p(z, E) we have 

inf d(x, zo) > llmzcritl + ^ (31) 

For simplicity we assumed that 9?7oo \ Vq just contains two complex conjugate 
points, but the extension to finitely many of them is straightforward in principle, though 
some of the estimates of Section S] would get more messy. Moreover, slightly more general 
perturbations g could be allowed, but we do not consider it worth the effort to include 
them. The conditions on /3 ensure that Zcrit is a critical point of p such that lim^^^^^j^ ^(z) 
exists. Finally, ((3T|) ensures that there is no sequence of singular points converging to the 
boundary of C/qo and is automatically fulfilled if e.g. p has only finitely many singular 
points. 

The integrability condition we need on p can also be expressed in terms of the metric 
d given in ([50]) . Let us define /i : R x A ^ [0, 00] by 



h{x,E) — sup 



dxPiy, E) 



p{y,Ey 



ye CR„+six) 



where A C M is the relevant set of energies, i.e. the support of the function Q{E, e). We 
assume: 



(A3) There exists S > such that 

/ h{x, E)p{x, E)dx < cx) 

jR\[a,fc] 

for some a < 6 G K. 

Proposition 1. Fix E <eR. Assume thatp{.,E) fulfills (Al). For x e E define £,{x) by 
iflD)) . and define p through p{£,{x, E), E) = p{x, E). Put 9'{^, E) = |§^. 

(i) Assume thatp{.,E) fulfils (A2). Then 9' fulfills Assumption 1 above on each com- 



pact subinterval ofM. containing £^^, with 7 — iJtt, 

^E) = Re(^crit), UE) - Im(ecrit) where Ccrit(^) := Uzcrit{E), E), 
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and for arbitrary a > 0. 

Moreover, 6' fulfils Assumption 2 on each compact subinterval of M such that 
inf{|x — ^rl '■ X <E 1} > 5, where 5 is the same as in (A2). In this case, k = \J S^ + (^ 
in Assumption 2. 

(ii) If p{., E) fulfils (A3) then ^ i-^ ll^'ll(f i £ -^s) ^^ integrable outside of a compact set. 

Proof, (i): We omit the dependence on E from the notation. By (A2) and (j28p . 

e(z)-Ccrit(z) = / {y-z,,itf{A + g{y-z,,it))dy = ——{z-z,,,t)^+\\ + 0{z-z,,it)). 
Moreover, 0'{O = ^ = ^, and thus 

Together, we get 

^'(^) = /i. \ (1 + ^(^ - ^-it)), (32) 

which gives Assumption 2. Now if (5 > is as in Assumption (A2), and if / is a compact 
interval with inf{|a; — £,^\ : x & 1} > 5, then p2|) together with Darboux's theorem for 
power series (cf. |BeTe2] ) yields Assumption 2 with k — yJQ + (5^ — ^c- 
(ii): By the results of [BeTe3| . it is enough to ensure that 

h{(,)=snv{\e\z)\:\z~-^\<i,+5} 

is integrable over the real axis outside of a compact interval. By the definition of ^ and 
the fact e'{i) = p'[t)/p^{t) this follows from (A3). D 

3 Time-Dependent Description of a Reflected Wave 

As explained in the introduction we construct solutions of the time dependent Schrodinger 
equation ([T]). Since we are interested in the scattering situation, we will need strong decay 
assumptions on the derivative of V , which we will give shortly. For the moment we will 
just assume that V is bounded on the real axis, analytic in a neighbourhood of the real 
axis, and that lim|3.|^oo dxV(x) = 0. 

A solution to ([T]) can now be given in the form of a wave packet: for any fixed energy 
E > maxa;gR y(a:;), let (j){x,E) be a solution to the stationary Schrodinger equation 

-e^dlc^=[E-V{x))c^ = p'{x,E)c^, (33) 

where 

p{x,E) = yjE - Y{x) > 

is the classical momentum at energy E. By standard results on ODE we can choose 
(/)(x, E) so that it is regular in both variables. Then the function 

^{x,t,e) := / Q{E,e)e-'''^'^(j){x,E) dE, (34) 
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for some regular enough energy density Q{E,e) supported on a compact set A, is a 
solution of H]). Of course, the solution to ((33)) is not unique unless we impose boundary 
conditions, and it is these boundary conditions that will determine the physics that the 
solution \1/ describes. In our case we will require that the wave packet is incoming from 
X = +00, as we will detail below. 

In order to connect to Section [21 we put 

and observe that (j> solves p3p if and only if ip solves 

ie9,V(^,i?)= (p2(f £;) l^^{x,E) = H{x,E)ijix,E). (35) 

Of course, this is just equation ^. 

Let us now fix our conditions on the potential V, in order to apply the results of 
Section [2] and to ultimately control the large x behaviour of the solution to ([55)1 . We 
select a energy window A = [Ei, £^2]. The most basic assumption is 
(VI) sup^ y{x) < El, and V is analytic in a neighbourhood U of the real axis. 

This assumption corresponds to (Al) from Section[2l For the condition corresponding to 

(A2), we define Uoo{E) as in ([29]) . We then assume 

(V2) Condition (A2) holds for Uqo(E) and for each E e A. 

Let us see what the somewhat implicit condition (V2) means for the potential V. Since 
p{x, E) > on the real axis, the singular points of p{z, E) will occur when V{z) = E, 
or when V itself has a pole. Condition (V2) states that for each E there is exactly one 
point closest to the real axis (in the metric d given in (|30|) ) where this happens. However, 
in contrast to the situation in the previous section, we now have to deal with a whole 
family 2;crit(£') of such points, indexed by the energy E. Therefore it might happen that 
for different values of E, different exponents f3 in (A2) occur, and E 1-^ Zcrit{E) needs 
not be continuous. In order to avoid such awkward situations, we impose 

(V3) The function E 1^ ZcritiE) is C"^ on A, and P does not depend on E in (V2). 

The most natural situation is when ^(^crit) — E with ^'(zcrit) 7^ 0. Then Zcrit(^) 
depends analytically on E. Another interesting situation is when V has a pole; then 
Zcrit{E) does not depend on E, but ^crit does via pI7)) . 

Finally, we need the analog to the integrability condition (A3), uniformly in E; indeed, 
in view of some L^ estimates to follow later on, we will need quite a bit more than (A3). 
We will say that V fulfils (V4)„j with ?7i > if the following holds: 

(V4)ni Define Ro{x) and Ciig{x) as we did just before ([29]) . We write 



x'^ 



{x) = Vl + 
and assume that 

h{x,5):= sup sup V'{z){E~V{z))-^/' < ^ (36) 



EeAzeCRo+s{x) 
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for some 5, z^ > 0, and all x outside a compact interval. Moreover, we assume that there 
exists C > such that for all sufficiently large a; £ M we have 

\V'{x)\ <C|x|-™-5/2 

Note that (V4)m implies, uniformly in i? G A, 

\p{x,E)~p{oD,E)\r^{xy^''^^"' as|a;|^oo. (37) 

Note also that since p(a;, E) is bounded above and away from zero uniformly in x and E, 
(A3) follows from ^. 

Given the conditions on V , we can now apply the results of Section [2l Recall 
the natural scale (^5]) and define p through p{^{x,E),E) = p{x,E) and 9'{^,E) = 
d$^p{£,TE)/p{S,,E). By Proposition [1] and Assumption (V2), 9' has exactly one pair of 
first order poles closest to the real line, which we denote by £,r{E) ± \S,c{E). Now it 
follows from Theorem [1] that in the n-th superadiabatic representation equation ([55]) 
reads 

/ Pn{e,i,E) £"+lfc„(£,e,^) \ 

with pn and A;„ as in Theorem [TJ Let us write 

Ci(£,X, S) 



C2{e,x,E) J -Me.i{x).E)^{T^4'){e,i{x),E) 

for the solution of ([55| with initial condition as in Theorem [31 If we choose n^ = ^c/e 
so that we are in the optimal superadiabatic representation, and translate back to the 
actual solution of ([55)1 . we find 

0(^,i^)=(7^„-/(';{^;^;;^j )) =ci(e,x,£;)^«+C2(e,x,£;)^«, (39) 

where the superscript means that we take the first component, and Lp^,j are the "optimal 
superadiabatic vectors" , i.e. the columns of T„^ . Note that due to ((23l) and (V4) we have 

^W(e,x,£;) = _i_ +0(£ < X >-(3/2+'')). (40) 

y/p{x, E) 

So by switching to the n-th superadiabatic representation we can write the solution of 
(|33p as a sum of two terms, which at the moment appear rather arbitrary. However, 
we will see now that, with the above initial conditions, the first one corresponds to a 
transmitted wave while the second one corresponds to an exponentially small refiected 
wave. 

To this end, note that by Theorem [3] we obtain 



Cl 



(e, X, E) = e^ ^» P^'-^^) ^'■(1 + 0(e)), 
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and 



^-a„-% 



^2^ \^ VMAE) J 

(41) 
where 7 = /3/(/3 + 1) is independent of i? by (V3). In the next section we will have to 
investigate x-behaviour of the error term very closely, but for the moment we do not 
worry about it. 

Since the vectors (p^j approach the adiabatic vectors as |x| — + 00, they are, to 
leading order, non-oscillatory at infinity. By considering the oscillations of the Cj and 
taking into account the fact that ^{x,E) approaches a linear function at x = ±cx), we 
find that the solution (|5^ corresponds to a wave incoming from the right at x ~ +00, 
described essentially by 

p-iI+°°ipiv.E.e)-pi+oo,E))/e 

^ ^-tp(+oo,E)x/e U2) 

\/p{'X),E) 
which is scattered on the potential. The refiected wave, for large values of x behaves as 

„iJ^'^ipiv,E,e)-p{+oo.E))/e 

\c2{e,+^,E)\- e»P(+oo,i=;)x/e_ (43) 

y/p(oo,E) 

When integrated against a suitable energy density Q{E,e) like the one we will specify 
below, this refiected wave gives rise to a freely propagating Gaussian, for large values of 
X and t, as can be shown by the methods of [JoMaj . 

Here, we are interested in the full history of the exponentially small reflected wave, 
which is only a part of the full solution that we would get by plugging (|5^ into (IM)) . 
From the discussion above we see that 

Xie,x,t) := / Q(£;,£)e-**^/"c2(e,x,^) ipi%,x,E)dE (44) 

gives the correct asymptotic behaviour of a reflected wave, and thus we define: 

Definition 2. The function x(x, t, e) as given in (|44)) is called the reflected part of the 
wave in the optimal superadiabatic representation. 

Let us now formulate the main result of this work. We start by fixing the assump- 
tions on the energy density Q. We assume that Q is given by 

Q{E,e) = e-^(^)/"e-*^(^)/^P(£;,e), (45) 

where: 

(CI) The real- valued function G > is in C^{A) and is independent of e. Moreover, the 
strictly positive function 

M{E)^G{E)+UE) (46) 

has a unique non-degenerate absolute minimum at E* g]£'i,_E2[. This implies that 
M{E) ^ M{E*) + ^'^"(^*) (^E - E*f + 0{E - E*f, with M"{E*)>0. 
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(C2) The real- valued function J is in C^{A). 

(C3) The complex- valued function P{E,e) is in C^(A) and satisfies lim^^o P{E,e) = 

P{E, 0) and 

"' "^ < C„, for n = 0, 1. (47) 



sup 

e>0 



dE' 



P{e,E) 



The first of our results shows how to approximate C2 in (|^1)) up to a small relative 
error. It should be viewed as an extension of Theorem [3] that includes control of L^(dx) 
norms. Let us introduce the shortcut 

y/2UE) V^UE) 

for the e-independent part of the argument of the error function in (|41l) . The point Xr{E) 
actually marks the intersection of the Stokes line emanating from £,c{E) with the real 
axis. Let E* be as in (CI), and let us write x* = Xr{E*). The point x* is the birth region 
of the reflected wave. We define for < J < 1/2 

( lix-x* < -ei/2-<5^ 

c,sis,x,E) = J 2sin(^)e-^4^e*«(-.^)-*«-(^)erf (^) ii \x - x;\ < e'/^-\ 
[ 2sin(^)e-^4^e*«(-'^)-i«'(^) ii x - x* > e'/^-\ 

and 



g{e,x,t) := f Q{E,e)c-''^/'c,s{e.x,E) ^l^l{e,x,E)dE. 
Ja 



Theorem 4. Assume that V fulfills (Vl)-(V4)m with m > 6, and that Q fulfils (Cl)- 
(C3). Then there exists C > 0, eo > and 5 > Q such that for all e < eo and all t € M 
we have 



|lXoff(e, .,t)~ x(£, ■,t)\\mdx) < C'e" 



e-i 



4+5 



The proof of this theorem will be given in the first half of Section U) Knowing now 
how much Xcff differs from %, we will need to show next that XcS itself has an L^ norm 
greater than that error term. Indeed, we are going to show more: we will give the precise 
asymptotic shape of the reflected wave for almost all times t after its birth time. To do 
so, let us introduce some additional abbreviations: 

R{E) -.^ -UE) - JiE), S{E,x,t):^-^ax,E)-RiE) + Et, (49) 

For an explicit description of Xcffj we have to consider three different regions of the 
variable x. The first formula will be valid near the birth region of the reflected wave. We 
define 

Xncm {£, X, t) ■- ^/ePa{x,E*)c ■= X 

fa{x,E*) \ _ 2 M"(E-) + iS"(E*.x.t) i=S'(B*,x,t)y , 

X / erf ^ — j^—^ + a [x,E )y \ e '' 2 e v? ^ '^ dy. 
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with 

Poix,E*) = 



2P{0,E*)sm{^) 



Here and henceforth, primes denote derivatives with respect to E. 

In the region of moderately large x, but away from the transition point, the reflected 
wave will be given by 

, , Po(x,E*)\/2TTe M(E-) + iS(i5*.x.t) S'(E-\^t)^ 

^M"{E*)+iS"{E*,x,t) 

For the asymptotic regime, i.e. for large x, the reflected wave is characterized by its 
asymptotic momentum 

k{E) := y^E-V{oo). 

rather than its energy. We put k* = k{E*), write /(fc) = f{E{k)) for any function /, and 
define 

u;iE)= / {p{y,E)-p{^,E))dy. 
Jo 

The reflected wave in the asymptotic region is then described by 

X{.,{s,x,t) = 2sin(7V2)P(0,g*)2V2^£fc* ^_MlEll ^^^inE')+u.(E')+.k'-E't) ., 



^M"{k*) + i{2t -{R + u)"{k*)) 



xexp i^-2kH+iR + ^nk*)r \ 

\ 2e{M"{k*)+i{2t-{R + Co)"{k*)))) 

Let us now give a precise statement about the shape of the reflected wave. 

Theorem 5. Assume that V fulfills (Vl)-(V4)m with m > 10, and that Q fulfils (Cl)- 
(C3). LetO <S < 1/2, Co > 0, and C\ > 0. For < l3 < 1/8 define 

{0 if X - x,{E*) < -Coe'^/'^-^ 

Xnc.Ae,x,t) tf\x-xAE*)\<Cos'/'-' 
Xfar(£,a;,f) ifCie''^<x-xAE*). 

Then there exist j3q < 1/8 and Eq > and D > such that for all (3 with (3o < (3 < 1/8 
and all < e < So we have 

llXcxpi(£, . , t) - Xcff(e, . , i)IL2(d,) < i^e3/4+.0e-«(^-)/^ 

for some Sq > 0. 

It is obvious that Xia.r is a Gaussian wave packet. Our next result shows that this 
is also true, to leading order, for all other times just after the transition occurred. The 
center of the wave packet is the classical trajectory qt of the reflected wave, where qt 
is defined to be the unique solution of the equation S'{E* ,qt,t) = 0. In particular, qt 
satisfies i J„^* , ^^^, dy + R'{E*) — t, which shows that qt behaves like t — R'{E*). 
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Theorem 6. Define 

Pq ( E* , qt) V2Tre M(g*)+iS(E'.x.f) ~ff,~l'^'^ „ 

' VM"{E*) + iS"{E*,qt,t) 

Assume t is such that e^'^^ < |gt — j:*| < e~^, with f3 < 1/8 as in Theorem\^ Then 
there exist C > 0, eo > and (5o > such that 

||Xgauss(£, . , i) - Xcff(£, . ,i)|li2 < Ce-^^ £3/4+^0 

/or aZZ e < £o • 

Like the two previous theorems, we will prove Theorem |6] in Section |4l A direct 
calculation of Gaussian integrals now yields 

Corollary 1. With the definitions above, we have 

4|P(0,£:*)sin(7r7/2)|7r3/4 m(e-) , 

||Xgauss(e, ■Ml^= IIXfar(£, ■ M L^ = \m"{R*)Y/^ "" ^^ ' 

Thus, for all times such that e^"^ < qt — x* , Theorem\^ yields an approximation Xcxpi 
of the reflected wave x whose norm is larger than the error term. 

Remark: We make use of the identity M"{k*) — AI"{E*)4:k* , which follows from the 
fact that E ^ M{E) is critical at E* . 

4 Rigorous analysis of the reflected wave 

In this section we give the proof of Theorems [3] - [SI It is checked directly that 

-oo 

(52) 
solves the second line of ((38|) with initial condition as in Theorem [3] In what follows, we 
will always understand n to be the optimal one given by (j20[) . 

The additional difficulty of the present problem when compared to [BeTell IEeTe2j 
is that we now need to extract the leading term of ([5^ with errors that are not only 
small pointwise in a;, but even integrable over all a; S M. In particular, neither the leading 
term of C2{e,x,E) nor all of its error terms are small when x — > +oo, and it is only 
after integration over E ^ A that these terms become integrable, due to their oscillatory 
nature. This determines our strategy: we will always decompose the error term into a 
part that decays as a; — > oo and one that oscillates there in a regular way, both of which 
will give rise to a small contribution to the i^-norm. More precisely, we introduce 



Definition 3. Let us write (x) = Vl + x^. 

(i) We say that a function r — r{e,x,E) is in £i{(3,i') for /3 G M and i^ > if there 
exists C independent of E such that 



|r(e,x,^)| <Ce'^{x)-'^ 



•2-v 
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(ii) We say that the function r is in £2 iP) for /3 e K if there exists C < 00 independent 
of E and ro(e, E) independent of x and with |ro(e, E)\ < Ce^ and 



, , , ei-P(-'^)ro(£,i?) if.T>0 

r[e,x,i^)-<^ if.T<0. 



Here, we put p(oo, i?) = Hm^.^ooPl-'c, S) = linij.^00 v^S — V{x). 

We immediately show how to use these two different functions in order to estimate 
L^ errors. 

Lemma 1. Let A be a compact interval and f : A —>■ M be a bounded function. For 
r = r(e,x,E) define 



I{x)^ / f{E)r{e,x,E)AE. 
Ja 

If r £ i?i(/3, v) for some u > Q then there exists C > such that 

mm^.)<CeP\\f\\,^^^E) 
If r Cz £2 (/?) then there exists C > such that 

ll^llL2(d2;) - ^^ \\f\\L^{dE) ■ 

Proof. The first assertion is just an appUcation of Fubini's theorem together with the 
fact that a function from £{(3, v) is square integrable in x. The second assertion follows 
from the fact that with 

J{x)^ I e-^P''°"'^y-'f{E)dE 
Ja 

we have 

f \J{x)\^ dx = ine f p{oo,E)\f{E)\^dE. (53) 

To see this, put k = p(oo, E) = ^yE -V{oo), so that E = P + V^(oo) and d^ = 2k dk. 
Putting A — p{oo, A) we have 



J{x)= / ei''='l^{k)f{k^ + Vioo))2kdk = V2^{T-^p){x), 

Js. 

where p{k) ~ 2fcl^(fc)/(fc^ + V^(cx))) and J> is the rescaled Fourier transform 

iT,g){k) = -^ f g(x)e-^^-l' dx. (54) 

By the Plancherel formula we find 



IK'pWUd.) = \\p\Wk) - 4y^ k'\f{k' + l/(oo))|2 dk = 2j^ ^E -V{oo)\ f{E)\^ dE. 
Taking into account the prefactor \/2'Ke we obtain (I53p . □ 
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Let us now start to see how to decompose quantities of interest into parts from £i 
and £2- 

Lemma 2. Assume that V fulfils (V4)m for some m > 0. Then for each /3 > and each 
< 1^ < m we can find ri € £i{P{m — v) ~ l,^) and r2 € ^2(0) such that 

l{a;>e-^}e*^^''^-^^ =ri+r2. 
In fact, we can choose 

with io{E) = J^ip{y,E)~p{c^,E))dy. 
Proof. We decompose 

^{x,E)^2 p{y,E)dy^2 {p{y, E) - p{(x^, E)) dy + 2K{x, E) + 2xp{(x,, E), 

Jo Jo 

where 

/•OO 

K{x,E)= / {p{y,E)~p{^,E))dy. 

J X 

By (V4)„i we have \K{x, E)\ < C(a::)^^/^^™ for x > 0, and from the inequahty |e'^ - 1| < 
\x\ valid for all a; € M we get 



1 — e2< 



■K{x,E) 



^ / \ — 1 /9 — ^ ^\^) 



-m+v 






Since (x) > e ^ by assumption, we have (a;) ™+'' < e'^^™ ''^ and thus 

(*) < C£'3(™-'')-i(a;)-i/2-'', which proves the claim. D 

Corollary 2. Assume that V fulfils (V4)m for some to > 0. Then for each fixed D ieR, 
and for each < ly < m, there exist 

^i^^ii-Trr-'A ^ 7-2 e £2(0) 

\ 2 + TO / 

with 

Proof. For /3 > and x with D < x < e^'^ we have 

^{D<x<e-l^ <^{D<x<e-P}{^} [x) <Ce ''"■' ' {x) ' 

which is in £i[~(3{l/2 + v),^). From Lemma [2] we know that the part with x > e^^ is 
in £i{(3{m — i^) — l,j/) + £2(0). Optimising over /3 gives the condition — /3(l/2 + jy) = 
/?(to - 1/) - 1, whence (3 = 1/(1/2 + to) and -/3(l/2 + i^) = -(2 + iy)/{2 + to). D 
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Corollary 3. Assume that V fulfils (V4)m for some m > 0, fix D E M. and < i' < m. 
Assume that f{e,x,E) is such that 



/oo 
f{e,x,E)dx 
-oo 



< Ce^, and Ix>d 



f{e,x,E)dx 



?(^ 



<C£^(x)~^-" (55) 



for some C > and some /3 G K. Then there exist 

2 + iy 



neSilp 



2 + m' 



^ , ^2 e £2{(3) 



such that 

Proof. Simply write /_^ f ds — J^^ / ds — /,°° s / ds, and use Corollary [21 
We will now apply the above tools to the study of ((52)) . Note that 



r(,(x.E) 

-«("'^) / f{e,s,E)ds^ri+r2. 

-OO 

pOO 

k(x 



lim 



lim yjE - V{x 



D 



(56) 



x^±oo X a;^±oo 

which is finite by assumption, so we can replace x with ^ for studying the decay rate. 

Proposition 2. Assume that V fulfils (Vl)''(V4)m with m > 0. For D sufficiently 
negative define 



a^) 



.z,(£,a:,i?) = l{,>c}-e*«(-) / e-iyie^+'k„^^ie,y,E))dy 



C2, 



For each z^q > there exists < ly < vq, and 



such that 



"^^^^'^-It^'")' '^^^^^Q 



\c2,Di£,x,E) -C2{e,x,E)\ ^ (ri + r2)e 



Proof. First note that by (V4)„i and ((22)) . for sufficiently negative y we have 

\e"+^kn{e,y,E)\ < e"(«=+«='"(«/«^»/^(2/)-(3/2+.) < e^^e"^^ (2;)-(3/2+'^). 
Using this in (|5^ gives 

|1{..<D}C2(£, 2^,^)1 <Ce°°e- 



-^4^M-(V2+.) = .-^4^, 



with r G i?i(oo, i^). Now let us define K{e, x, E) = J^ {Pn{y) — 1/2) dy which is finite and 
0{e^) for all a; by ([21]). We write 



?(:>=) 



P«(y)dy=-C(a;,i?)-i^(£,e(.T),i?)+if(e,0,i?). 



(57) 
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Now for X > D we have 

/•OO 

\K{s, ax),E)\ < s' / (y)-3/2-- dy < C£2(x)-l/2-^ (58) 

by (HH). Thus 

l{.>D}e~^^(^'«(^)'^) = 1{.>D}(1 + £r(:r, £, E)) (59) 

with |r(x,e,£;)| < {x)-^-". Using dH]), ([21 and (V4), we find 

|e"+iA:„(e,:E,S)| < CVie-^(.T)-3/2--^ (60) 

and ([5^ together with ([SO]) give 



— oo 



with ri e £i(^, i^). Inside the integral, (P7)) gives 

e-i /_"^ ''"W'i>+ ,„(y, r) = e"^ ^— p^Mdr^^ ^ 0(e°°)), 



and using ((2T|) we see 

e-f /-"oo P.M dr+iK(e,o,B) = e-ii. + 0(g)_ 

Combining this with ((60| and ((57l) we obtain 

ei(/-"» p™('^) <i'-+^(-o.^))e"+iI„(e, y, i?)V+,4£, V, E) = 
= e-i%"+ifc„(e, y, ii;) + £3/2g-^^(^^ ^^^ ^)^ 

where |/(e, j/ji?)! < (a;)~'^/^~*. When inserting this into ([61]l . the first term above gives 
C2,n, while the second one is estimated using Corollary [3] to give the desired result. D 

In order to proceed further, we now need the following refinement of Theorem [2] 
which gives us additional control over the x-dependence of the error terms to fc„^ . 

Proposition 3. For each a > there exists r e £i(3/2 — a^l + v) such that 

£"^+ifc„, (e, e, E) = g{s, e, E) + e~^r{e, C, E), 
where g is given through 



Proof. We will use Assumption (V4). We use ((25|) in (finite) neighbourhood of ^r and 
([m outside this neighbourhood; As for ([^^ . we know ([BeTeSj, Proposition 2) 

and 9'{^) — p'{£,)/p{£,)i together with (|56p and (V4),„ the claim follows. Note that we 
may take a arbitrarily small in Assumption 2, due to Proposition [TJ □ 
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We will also need the following useful little lemma: 
Lemma 3. For each n G N and each e > Q we have 

a;"e ^^ g £i{n/2,v) for all v > Q, 
and, for (5 > 0, 

l{2;>£-*}6 ^^ G Si{m, v) for all m > 0,1^ > 0. 

2 

Proof. Taking the derivative of g{x) — x"e~T^ and equalling it to 0, we see that g is 
maximal at x = \j2en. The value at the maximum is given by e~"'^(2en)"'^ = ©(e"'^). 
On the other hand, f{x) = e~4F < e~^ < C{x)~'^ for all i/ > 0. Since a;"e~^" = 

,2 

f{x)g{x) and e~ 27 := f'^{x), both statements are now obvious. D 

Proposition 4. Assume V fulfils (Vl)-(V4)m- Define 



C2 



Lie, X, E) = l{,>c}2 sin (^) e-^e*«-)-t«'(^)erf 



T/ien for arbitrary a > 0, for D negative and \D\ large enough, there exist 

such that 

C2{e,x,E) = C2xie,x,E) +e^~^{ri +r2). 

Proof. By Proposition [2] it is enough to consider 02, d — C2,l- Using Proposition [3] and, as 
before, ([60|) and Corollary [3l we arrive at 

C2,z) = l{x>i5}-e*«("W e-iyg{e,y,E)dy + c-¥{n+r2) (62) 

^ J — 00 



and ri, r2 as in the statement. Now recall that 

ff(£, y, E) = -2i JJ sin (f^) e"^ e"^"^ cos {w{e, y, £;)) , (63) 



with 

( T^\ y^^' Jy-^rf , (r,{y-Q 

-^('^y^E) = — ^^^ + ^^^- 

We write 

cos {w{e, y, E)) = - U'^i'^^y^E) _^ ^-iw(e,y,E)\ (g4>) 

For the second term above, the exponential e~^ present in e^™^^^^'^^ combines with the 
one present in (|62p : we use integration by parts in (j62p . with the function e~^^ being 
integrated. The result is 



e '''e 



Ly i^(e,y,E)-^^i^^ ,„ _ ^^-i?^-™(e,C,-E) -^^^^ 



e — 2^^- dy = — :e-?^e-™^'''^'^^e — 5^«^- (65) 



-2i 

^ f e-i^e-™(---)e-^^ f^i^^ill - ^^^^ + i.^ 

2l J-00 \ Sc $c s< 
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The boundary term is clearly in £i{l,iy), and the factors left out in ([^5)1 are bounded by 
Ce~^/^e~^, so together this gives a contribution in e~^fi(l/2, i^). The contributions 
of the integrand in the bulk term stemming from the last two terms in the bracket are 
in £1(1, 1/) by Lemma [31 For the first term, we have to integrate by parts once again; the 
resulting boundary term is again clearly in £i(l, i^), while the bulk term is very similar 
to the one in (j65p with the crucial difference that the bracket there is now multiplied by 
2(y — ^i.)/'fc and an additional term 2e/^c is present. Applying Lemma[3]then shows that 
this integrand is in fi(l, v) again. Now we can apply Corollary [3] and find that the term 
involving e^^^^V'^'' is bounded by e-^(ri + rg) with n e £i(l/2, ly) and rs G ^2(1/2). 
We now turn to the term involving e+™'^'^'^'^^. Here the oscillations e~^ cancel in the 
leading term of (|62p . which is therefore given by 






2,D = h.>D}Jj^sm{^)e~^e^'^^^>-^<' e^^^e ^^ ' ^" dy. 



i(x) 



iv-ir)^ _i(«-er)%i'^e(y-S,) 



We now use explicit integration 



/ e-^J«^d2/=V^^erf^-=4L 



and the fact that |e'^ — 1| < |a;| for all real x to find 

C2D =c2,L + l{x>D}^e - e' - / e 2.5, f{y)dy, (66) 

with 

Using integration by parts, we find 

/'' e-^ ^ dy = 2e£e(2 - e-«'/(2^«=)) - g^^VC^eC.) j^2 
■/-oo eCc 

Now an easy calculation using Lemma [3] shows that Corollary [3] applies to the second 
term of (pS)) . showing that it is bounded by e~~^{ri + r2) with ri e £i(l/2,j/) and 
7-2 e £2(1/2). D 

After all these preparations, we now work towards the proof of Theorem U) Recall 
that X was defined in (j44|) . and in a similar way define 



XL{s,x,t) := / g(^,£)e-^*^/^C2.L(e,x,£;) ip''^'l{e,x, E)dE. (67) 

Proposition 5. Assume that V fulfills (Vl)-(V4)m, and that Q fulfils (C1)-(C3). Then 
for all a > there exist £0 > and C > such that for all e < Eq we have 

||x(£, . , t) - XL{e, ■,t)\\^.< Ce-^^£i-^-". 

In particular, for m > 6 we have ||x(e, . ,t) — xl{£, ■ , OIIl^ < Ce s £3/4+5 j^^ some 
S>0. 
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Proof. By Proposition |3] and Lemma [1] we have 

1^2+^-r. , (1), 



||x(e, . , i) - XL{e, .,t)\\L^< Ce-^~^-" sup |v'i,il(IIQollLi(d£) + e'^' HQoWlh^e))^ 



x,E 



where Qq{E) = Q{E)e "^ . (|^ ensures that sup^,^ |(^^ 2(2: 2;, i5)| is bounded uni- 
formly for aU e < £o- By virtue of (CI), for arbitrarily small (5 > we have M{E)/e > 
M{E*)/e+e-^^ when \E-E*\ > e^'"^-^ . Thus writing A^ = {£; e A : \E-E*\ < £i/2-5} 
we find that for each m G N there exists C > with 



|Qo(i?)l{B^Aj|<Ce ^e-, 



(68) 



and, by (CI) again 



|Qo(i?)l{£eAj| < Ce ^e ^ (i+(^ ». 



Consequently, by scaling 

T^riT^*\l f m" (E*)iE-E*)^ -, M(E*) 

IIQoIIli < Ce^'"'^^^'^ e -^^ -^dE + C\^\e V^e™ 

< Ce-*'^(^*)/^ei/2. 
Similarly, ||Qo|Il2 < Ce-'^'^^'^'^l^e^'^. This shows the claim. 

Formulas (|68|) and ((69)) are also the key to the following 
Proposition 6. Assume that V fulfills (Vl)-(V4)m, and that Q fulfils (C1)-(C3). Put 

Ae = {£;e A: |£;-£;*| ke^'"^-^}, 

and 



(69) 

D 



„(i)/ 



XL,o(e,x,t) := / Q{E,e)e-''^^'c2,L{£,x,E)^\'>,{e,x,E)dE. (70) 

J \^ 

Then for all m > there exist Eq > and C > such that for all e < Eq we have 



|XLo(£,-,t)-XL(£,.,t)|lL2 < Cc 



M{E*) 



E . 



Proof Note that by (gO]) we have (^^^^(e, a;, £■) 
moreover, we have 



1 



00 



ri with ri e fi(l, 1 + i'); 



1 - erf (u;) = ^ / e"" du < ^ / -e"" du = ^^e-""" (71) 

for w > and erf(u') < 26""" /(Trjuil) for w < 0. This shows C2,L(e,x, i?) 1^9^ 2(272;, i?) = 
e~"J^(ri + r2) with ri G fi(0, 1) and r2 G ^2(0). Then (|55|) and Lemma [T] shows that 

g(^, e)e-"^/=C2,L(e, x, S) ^f ^ (e, x, S)d£; 



A\A, 



< Ce V^e" 



Li 



for any m, the same holds for the corresponding L^ norm, and the claim follows. D 
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The final step in tlie proof of Tiieorem [5] is an effective description of the error 
function in C2,l. Define Xr{E) and a(x, E) as in (|48)) and note that E i-^ x^{E) is C^ by 
(V3). Put 

■/ p\ - / 1 if a; > x^{E), 
J{x,i^) I Q [ix<x,{E). 

Proposition 7. For any m > 0, fc > anrf (5 > we /laue 

sup \eri(a{x,E)/^/e) ~j{x,E*)\ In^-^jE'Mye^/^-^n € £i{m,k). 

EGAe 

Proof. Let us consider only the case a; > x^{E*), the other case goes analogously. We 
first show the estimate x — Xr{E) — {x — XriE*)){l + 0{e^)). By (V3) we have \xi-{E) — 
Xr{E*)\ < C\E - E*\ for some C > 0. Thus, 

X - Xr{E) = {x- x,{E*) - C\E* - E\) >\x- x,{E*)\ - Ce^/^-a 

on Ae, for another constant C > 0. Since |a; — x,:{E*)\ > ei/2-2<5 impHes e^^'^^^ < 
\x-x,{E*)\e^, we find 

X - Xr{E) = {x- x,{E*)){l - Ce^). 

By ^ we find 

j(x, E*) - CTi{a{x, E)/^e) = 1 - erf (a(a;, E)/^e) < ^ „^!(^^ 

Now since p{y,E) > c for all y and all F, and since S,ciE) is bounded away from 0, we 
get 

1 r 

a{x,E) = / p(y,F)d2/> co|a;-Xr(F)|. 

From the estimate above we find 

a{x, E) > co(x - x,{E*)){l ~ Ce^) 

on A^ X {a: : |x — Xr(£'*)| > e^/^^^''}. Plugging this into (TfSjl and using Lemma[3l we 
obtain the result. D 

Now Theorem 2] follows by using Proposition [7| in the expression for C2,l in Propo- 
sition[B]in case jx — a;*| > e^^''^~^ , while keeping C2,l for |a; — a;*| < s^/'^~^ , and then using 
Lemma [T] again. 

We now turn to the proof of Theorem [5l Let us recall the abbreviations 

M{E) = UE) + G{E), S{x,E,t) = ^E) + J[E) - i^(a;,i?) + M, 

and 

a(x E) - ^(^^E)-^riE) ^ Jl^E^piy,E)dy 
V^UE) V^UE) 
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We first take advantage of Proposition [SI i.e. the fact that we know E — E* is small in 
the relevant area of A, and expand all quantities of E around E* . The only exception 
is ^(x, E) since x i-^ ^{x, E) grows like p(oo, E)x, and so in this case such an expansion 
would not be uniform in x. We have 

±:^^±^^d,aixE*) nE-E*fX^ (73) 

\/e y/e ^Je \ y/e J 

where the error term is not uniform in x. By the Taylor formula with integral remainder 

erf(z + ^o) = erf(0) + ^ / e-(^+"^«)' ds, (74) 

Vi" Jo 

we can write 

„,(^)^„,(^^^,^,,_,.,^^„(,^^ ,,,, 

The error term is uniform in x here, because the error functions are. Since we are on A^, 
by restricting attention to the set {|a; — x*| < e^/'^~^} we actually get 

O (^^-^\ l{|._,.|<,v.-.} e £M'^-^\^)- (76) 



Hence, by (|69|) and Lemma [T] again, we find that once multiplied by Qq{E) and the , the 

L {dx)-norra of the corresponding dii'- integral of this remainder is of order 0(6 J e 
Similarly, we expand 

So{E,t) ■.^+£,r{E) + J{E)+Et 

and M{E) to third order. The error stemming from this expansion is of order \E— E* p/e 
and is uniform in t and x. Thus, by the same argument, this remainder eventually yield 
an error of order 0{e ^ e^~'^^) in L^((ia;)-norm on the set {\x — x*| < e^/^~*}. 

Finally, we use (|40)l to see that (p \{e, x, E) = . ^ + ri with ri e £i(l. 1 + i^) 

y/p{x,E) 

and estimate the error term using the by now familiar reasoning. Defining 

-,,-■ ^mi^) sin (.7/2) 

^ypix, E) 

and expanding around E* we obtain 

Po{x, E) = Po{x, E*) + 0{E - E*), 
with an error term that is uniform in x and of order e^"^ , so that we have shown 
Proposition 8. Define 

XL,i(£,x,i) =e-(*^(^*)+'^«(^*'-'*))/^/^ Fo(2:,^*)e-'^^T^^o(^''*)x 

(77) 
xe-^^^?^(^"(^*)+'^o'(^*'*))ei«(-^^)erf [^^^ + ^^^^{E - E*)) dE, 
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where all the primes denote E -derivatives. Then, for all < S < 1/4, 

To get the result stated in Theorem[5]on the range {|a; — a;*| < e^^^~^}, we need to 
expand the term ^(x, E) around E* as 

(78) 
where primes denote ^^-derivatives. Over the considered set of x's, the error term gener- 
ated is dominated by \E — E*\^ /e, and is deah with as above. To arrive at the expression 
Xneor, we change variables from E to y = [E — E*)/y^ in the remaining integral, and 
note that extending the range of integration to the whole of R costs an error of order 

gOOg-M(B*)/e^ 

We now consider the region e^^^"^ < \x — x*\ < e~^ with (3 < 1/8 and < S < 1/4. 
The following proposition settles Theorem [5] for the intermediate region. 

Proposition 9. Define Xmod cis in i5U\). Then there exist Eq > and C > such that 
for each e < Eq we have 



(Xmod - XL,l)l{£i/2-«<a;_a;.<e-^} 



L2(da;) 



< Ce-^^e^-^^-^^. 



Proof. The leading term is obtained from the leading term of (j78p : we know from Propo- 
sition [7] that we can replace the error function in (j77[) with 1 for the range of x we 
consider, and then we obtain a purely Gaussian integral in {E — E*). Again, we can 
extend the range of integration to all or R for the cost of an error of 0{e~^^^ ^/^e°°). 
When integrated over e^/^"'' < x — x* < e~^ this error stays of that order, and distribut- 
ing the various leading terms of ([75)1 in the exponent changes the 5*0 and its derivatives 
to S. For the error term we note that \E — E*\ < e^l"^^^ on A^, and thus 



for I a; I < e ^ . Since the length of the range of integration in a; is 0(e '^), the length 

of integration range in _E is e ' , and the factor e ^ is present, we obtain the 
result. D 

It remains to treat the region x > e~^ . Since this has essentially been done before 
I JoMaj ■ we will be much briefer here. We use Lemma [2] in order to write 

l{.>e-^}e*«^"'^^ = ei("(^)+P(°°^^)-)(l + n), 

with ri g £i{[3{m — u) — l,i^). By Proposition [71 we can replace the error function 
by 1, and replace p{x,E) by p{oo,E) with an error bounded by C/{x)™^^/^ and thus 
by C'e''(i+"-'^)(a;)^+i/2 « Ce^^"'~'''>~^{xy+^/'^. Integrating these error terms over E 
against the Gaussian function gives a term whose L^{dx) norm is bounded by means 
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of Lemma [TJ The L^{dE) norm of the Gaussian in energy is of order e^/^, see ([M]) . 
so that the error term is eventually of L^{dx) norm of order 0(e ; £;/3('"-f)-i/2')_ 
We need /3(m - v) - 1/2 > 3/4, so m > 5/(4/3) > 10 with (3 < 1/8. Finally, the 
leading term can be computed by Gaussian integration. First we change variables from 
E h^ p{oo, E) = ^E-V{oo) tok^ E{k) = fc^ + y(oo) and write f{k) = f{E{k)) for a 
function /, just as we did in the proof of Lemma[21 Carrying out the resulting Gaussian 
integration, we find that if m > 10, then there exists 5 > Q such that 

XL,0(e, X, t) = Xfar(£, X, t) + 0(e ^£3/4+5)^ 

where Xfar is given by ([5T|) . This finishes the proof of Theorem [5l 

It remains to prove Theorem [51 The idea of the proof is to observe that Xnoar and 
Xmod are sharply peaked around the classical trajectory qt of the reflected wave, i.e. the 
solution unique solution of S'{E*,qt,t) = 0. This will allow us to expand the reflected 
wave around its maximum in the variable x, and get the result. Since we already know that 
Xfar is a Gaussian centred around the asmptotically freely moving classical reflected wave, 
I JoMaj ■ it will be negligible for the times t considered. We start by proving concentration 
in X for flnite times. Recall that by the change of variables y = {E — E*)/ ^/e, we can 
write, for any Z? e M 

Af (g*)+ig(-E* ,:i:,t) 

XL,i(£>2;,t) = l{^>c}V^-Po(e,2;,£'*)e ^ x (79) 

X / erf +a {x,E )y \ e y 2 e v^ ^ ' '" dy. 

J\y\<e-^ V V^ / 



Proposition 10. For all m > 0, all S > and all C > there exists C > such that 



<Ce — 



for all t e R. 

Proof. Let us flrst note that we can extend the range of integration to all y G M in ([79]) at 

the expense of an error bounded by Ce ^ e™ in L^{dx) sense, by the usual arguments. 
For (5 > choose A; G N such that k > m/6. fc-fold integration by parts gives 



XL,i{£,x,t) = l{^>£,}VeFo(e,a;,£'*)(-l) e - x 

dy* 



^ ("erf h^^^2^ + a'{x,E*)y\ ^-y\M"iEn+^s"iE' .,.,t)),2\ ^ 
{-iS'[E*,x,t)Y (ly + uye e j. 



The fc-fold derivative has the form f(e,x,y)exp{—y'^M"{E*)/2), where 

\fie,x,y)\<Cik)\y\'' (80) 



28 



uniformly in £ > 0, in time t and |x| < C. To see this, note that by the formula dl.g{a + 
bx) — b"-g''''\a + 6x), the l-th derivative of the error function is an Hermite function, i.e. 
is given by 

a'(x, E*yHi^i {a{x, E*)l^e + a'(x, E*)y) e-('^(-^^')/Vi+a'(x,B*)y)^^ 

where Hi is the l-\h Hermite polynomial. Since Hermite functions are uniformly bounded 
in their variable, the uniformity claim in e follows. Since the derivatives on the second 
factor also produce Hermite functions, we have shown ((80|) . Moreover, \S'{E* ,x,t)\ > 
c\qt — x\ for some c > uniform in i, and so 

fe/2 



^{\i^-\>-''^-''^\S'{E\x,t)\ 



< e^yc" < e'^/c" 



Thus the dy integral is bounded by C{S, m)e"^ and bounded uniformly in x and t, and 
estimating the x integral by the area |a;| < C times its maximum gives the result. D 

In the region where x is moderately large but not finite as £ —> 0, we need a different 
argument due to the broadening of the wave packet. 

Proposition 11. Put 



A = {a; e R : e^/^'^ ^yM"{E*)'^ + 2S"{E* ,qt,t)^ < \x - qt\ < £-'^}. 
Then for each m d N we have 

M{E*) 

||XmodlA|lL2(d:r) < -ft^e ^ £", uniformly int. 



Proof. First recall that S"{E* ,x,t) is independent of t and grows with x like x. The 
negative of the real part of the exponent of Xmod is given by 

S'{E*,x,trM"{E*) 



2e{M"{E*)'^ + S"{E*,x,t)'^)' 

Since \S'{E*,x,t)\ > c\qt - x\ for some c> and \S'' {E* , x) - S'' {E\ qt)\ < C\qt - x\ 
for some C > 0, we find 

/(,) > C^\x-q,\Hl"iE*) 

■'^ ' - 2e{M"iE*)^ + 2S"iE*,qt,t)^ + 2C\qt~x\^)' ^ ' 

Differentiation of the function y i-^ ^^^ ^^ shows that the right hand side of ([5T|) is 
minimal when |a; — q^p is minimal, i.e. when |a; — gt| is equal to its lower bound in A. 
Inserting this yields 

.,. -. .-28 cHM"{E*f + 2S"{E\q,,tf)M"{E*) ^_,, ^M"{E*) 



2(Af"(£;*)2 + 2S"{E*,qu tf){l + 2C£i-25) 2(1 + 2C£i-2'5) ' 

A-/(£;*) 

This shows that |xmod| < Ce e £™ for each to e N all a; G A, with C uniform in t. 
Now, since for t large g* ~ i so that S"{E* ,qt,t) ~ t, we get that on A, t < ci£~^~^/2+'5^ 
so that \x\ < C2e~^~^^'^^^ , for some constants ci,C2. Integrating over the range \x\ < 
C2S^^^^I'^^^ then shows the claim. D 
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We can now complete the proof of Theorem [HI Define Xgauss as in Theorem [HI Put 
K^ = {x:x~x*< iei/2-5}^ ^^^ K^ = {x : x - x* > 2e-'^} 



llXgauss -Xcffll < ||Xgauss(lA+ + 1a_)|| + || Xnoar 1a_ || +||xtarlA+|| + 

(Xgauss - Xmod)l{£i/2-a/2<a;-a;;r<26-''} +0(^6" « £ 



0(e-^^^.3/4+.o 



where all norms are L^ norms, and x* — Xr{E*). Now for e^^'^ ^ < \qt — x*\ < e ^, 

the first and third term are O ( e « £"M for any tti > 0, as can be seen by direct 

estimates on Gaussian tails. By Proposition [TU] the second term is of the same order for 
\qt — a;*| > e'^/'^^^ . For the term in the second line, we may restrict our attention to the 
region \x - qt\ < e^'^-^ ^ M" {E*Y + 2S" {E\qt,tY . 

This follows from Proposition 1 11 1 for Xmod, and from direct Gaussian estimates for 
Xgauss- For these x we expand the exponent of Xmod! we have to investigate 



d{x, t) :- 



S'(E\x,tf 



{x - qtf 



s{M"{E*)+[S"iE*,x,t)) Aep{E\qtY{M"{E*) +iS"{E* ,qt,t)) 
Since both \dxS' {E* ,x,t)\ < C and \dxS"{E* ,x,t)\ < C uniformly in x and t, we have 

.2 {x~qtf 



S'{E\x,tY 



4pHE*,qt) 



<C\x-qt\' 



and 



{M"{E*)+iS"{E*,x,t)) iM"{E*)+iS"iE*,qt,t)) 

\x-qt\ 



<C- 



M"{E*)+iS"{E*,x,t)\\M"{E*)+iS"{E*,qt,t)\' 



Then we plug this into the expression for d(x, t), and make use of the inequality |a6— cdj < 
\a — c\\b\ + \b — d\\c\, and of the estimate 



M"{E*)+iS"{E*,x,t)^ 

= iM"{E*)+iS"iE*,qt,t)){l + 



{M"{E*) + iS"{E*,qt,t))[ 1 + 



i{ S"{E*,x,t)-S"iE*,qt,t)) 
M"(E*) + iS"{E*,qt,t) 
\x-qt\ 



\M"{E*)+iS"iE*,qt,t)\ 



(82) 



On the set \x - qt\ < V2e^/'^-^\M"{E*) +iS"{E* ,qt,t)\, the last term of the parenthesis 
of order e^"~ , uniformly in x,t. Eventually, we get on the set 

\x~qt\< e^/^-^^M"{E*f + 2S"{E\qut), 

for uniform constants Ci, C2, 



d{x,t) <— |x-gt|3<C2e 

£ 



1-35 
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It then remains to make use of these estimates in the computation of the modulus of 
Xgauss - Xmod On the Set \x - qt\ < e^l^-^ ^ M"(E*Y + 2S"{E* ,qt,t). We finally get, 
with a constant C3 uniform in i, 



l-3<5 



KXgauss — Xmod)l{£i/2-«/2<2:_:E*<2e-/3}| < C^\Xga.uss\s 

Explicit integration of the right hand side above finishes the proof. 

Let us finally remark that in a similar way, one can expand the expression (j79p for 
the near region, and get an approximate expression for the transition region as well. 

5 Superadiabatic representations 

In this section we will prove Theorem [TJ We start with equation (jlip which we repeat 
here for convenience: 

ied^^{S,) ^ HiOi^iO :- \ (^ J^^^^ ^ ) V(0- (83) 

While the ultimate aim is to construct the superadiabatic transformations T„ from the- 
orem [U we start by constructing superadiabatic projections. For this matter, write 

--(^)-Kp(O t)' -(0^1 -MO. 

for the spectral projections corresponding to the eigenvaluse ±1/2 of H. The n-th su- 
peradiabatic projection will be a matrix 



^(»)=^efc^,, (84) 



fc=0 



with ttq{x) ~ T^+ix) and 



(7r("))2-7r(") =C'(e"+i), (85) 

ie^e - iJ, 7r(")j =C'(e"+i). (86) 

Here, [A, B] — AB — BA denotes the commutator two operators A and B. Likewise, we 
will later use [A, B]j^ = AB + BA to denote the anti-commutator of A and B. Equation 
(|55|) says that tt^"' is a projection up to errors of order 6"+^, while ([5S)) states that if 
^ 1-^ i/;(^) e C^ solves ((83|) . then also ^ 1— > tt'-"' (0-0(0 solves ([83)) up to errors of order 

We construct tt^") inductively starting from the Ansatz ([M]) . Obviously ([55]) and 
((55|) are fulfilled for n = 0. In order to construct 7r„ for n > 0, let us write Gn{t) for the 
term of order £"+^ in ([55)) . i.e. 

(7r("))2 - ^(") = e"+iG„+i + 0(e"+2) , (87) 

with 

n 
Gn+1 = /.^T^jT^n+l-j- (88) 
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The following proposition is taken from [BeTel| and the proof there applies literally 
also to the present case. 

Proposition 12. Assume that tt*^") given by |g^[ ) fulfills i85\) and i86]) . Then a unique 
matrix tTu+i exists such that 7r("+^' defined as in ^84^ fulfills i85\) and i86]) . 7r„+i is given 
by 

TTn+l = Gn+1 - [Gn+i, TTq]^ - i [tT^, TTo] . (89) 

Furthermore tt^ is off-diagonal with respect to ttq, i.e. 

ttott'^tto = (1 - ^o)<(l - ^o) = 0, (90) 

and G„ is diagonal with respect to ttq, i.e. 

7roG„+i(l - vTo) = (1 - Tro)G„+i7ro = 0. (91) 

As in |BeTel| we represent the matrix recurrence relation (j89p with respect to a 
special basis of C^^^, in order to arrive at a more tractable set of scalar recurrence rela- 
tions for the coefficients. While the following matrices are different from the one chosen 
in [BeTel] . they satisfy algebraic relations very close to those obtained and utilized in 
|BeTelj . As a consequence most of our proofs require only slight adjustments as compared 
to those in [BeTelilBel^ . Let 

X^(-; J), Y = -2H, Z = :±Y', 

where 9' — p' /p. Together with the identity matrix 1, X, Y{^) and Z{^) from a basis of 
jj2x2 £qj, every ^ G M. The following relations now follow without effort: 

X' ^0, Y' ^ -e'Z, Z' = -9'Y, (92) 

[X, Y]^ = [X, Z]^ = [Y, Z]^ = 0, X^^Y^ = -Z^ = 1, (93) 

[X,Tro]^Z, [y,^o]-0, [Z,no]^X, (94) 

l-[l,7ro]+ = r. (95) 



The relations (|92|) - ((95)) agree almost exactly with the corresponding relations derived 
and used in |BeTelj . The only differences are a sign change in the last equation in ([^^ 
and in the last equation in ([M)l . These relations show that this basis behaves extremely 
well under the operations involved in the recursion (|89p . This enables us to obtain 



Proposition 13. For all rt G N, 7r„ is of the form 

7r„ = XnX + UnY + ZnZ, (96) 

where the functions a;„,y„ and z„ solve the scalar recursion 



Xn+l 
Vn+l 



iO'vn (97) 



y^ {XjXn+l-j + yjVn+l-j - ZjZn+l-j) (98) 

J = l 

Zn+1 = -ix'„, (99) 
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and satisfy the differential equations 

x'„ = izn+i, (100) 

y'n = ()'^n, (101) 

z; = ixn+1 +0'yn. (102) 

Moreover, 

MO = ^o'iO, yi(0 = ^i(0 = o. (103) 

Remark 1. Hence, for all even n, Xn = 0, while for all odd n, 2/„ = z„ = 0. Moreover, 
Xn is always purely imaginary, while j/„ and z„ are real. 

Proof. For n = we have Gi = and thus 

which shows (|103p . For the rest we just have to stick (|M|) into the recurrence ((5^)) and 
use the relations (P^ -" ([55)l . For Gn+i we use ([M|) to obtain 

n 

G„+i = ^{xjX + yjY + ZjZ){xn+i-jX + yn+i-jY + Zn+i-jZ) 

= X! i^3^n+l-j + y-jUn+l-] - ZjZn+l-j) 1 . 
J = l 

With ([M]) this yields 

n 

Gn+1 — [G„+i, 7ro]_|_ = ^^ (a;ja;„+i_j + yjyn+i-j — ZjZn+i^j) Y . 



Next we compute with (|92|) that 

< = x'„X + y'^Y + z'^Z - e'yr^Z - 0' z^Y , (104) 

which implies with (|94p that 

K,7ro]^x'^Z + z'^X-e'ynX. 
In summary we obtain 

n 

TTn+i = ^ {xjXn+i-j + yjy„+i-j - z^z^+i-j) Y ~ i [x'^Z + z[^X - 9'ynX) , 

from which dnT])-® can be read off. While PUO]) and (fTU^ follow directly from ^^ 
and dMl), for (fTUTl) we observe that 

= TToTrJjTTo = {-y'n + 0'Zn)TTo , 



where the first equality is (j90|) and for the second equality we use ttoXtto = ttoZtto = 
(a consequence of ([M)) ) and ttoFtto = — ttq (a consequence of dM])) in p04p . D 
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Remark 2. From (jlOOp through (|102p we may derive recursions for calculating x„ or 
z„, i.e. 

xn+2{0^-<iO + o'iO I^Je'iOx'^iOd^ + c^ . (105) 

and 

^n+2(0 --^ (4(0 -^'(0 (/e'(0^«(e)de + c)) . (106) 

The constant of integration C must (and in some cases can) be determined by comparison 
with (P7)) ~ ([M)) . In the case where 9' is given as in Assumption 2, this strategy will lead 
to fairly explicit expressions of the leading order of the coefficient functions a;„, ?/„ and 
z„, cf. Proposition [TBI from these we will extract the asymptotic behaviour of a;„, z/n and 

Using (|100p - (jl02p . we can give very simple expressions for the quantities appearing 
in ([55)1 and ([M]) . As for ([5S)) . we use ([5^ and the differential equations to find 

ied^ - i/,^(")] = ie"+V; = -e"+\zn+iX + Xn+iZ). (107) 

The term by which tt^") fails to be a projector is controlled in the following proposition. 
Proposition 14. For each n e N, there exist functions gn+i,k, k <n, with 

((^(n))2 „ ^(„))(^) ^ fJ2 S^+'^gn+lAOj 1- (108) 

For each k < n, 

ffn+i,fc = 2i(a;fcZ„+i - z^Xn+i). 

Proof. Recalling ([M)) we have that 

n n— fc 

fc=l j=0 

As before we compute 

n—k n~k 

y^ 7rfe+j7r„_j = ^ (xfe+j a: + Vk+jY + Zk+jZ) {x^-jX + jjn-jY + z-a-jZ) 

3=0 j=0 

n — k 
~ 2^ i^k+jXn-j + Uk+jUn-j — Zk+jZn-j) 1 =: ffn+l.fel ■ 
j=0 

With Proposition [T3l 

n — k 
9n+l,k = X! K^k+j+l^n-j + Xk+jZn-j + l) + {O' Zk+jVn-j + O'Vk+jZn-j) + 
3=0 
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— i{xk^j^iZn-j -\- Zk+jXn-j + l) — {0 yk+jZ^-j + Zk^jljn-j) — 
n—k 
— 1 / ^ \^k~\-j''rl^7i—j ^k-\-j'^n—j-\-l i '^k-\-j ^n—j-\-l '^k-\-j-^l^n—j } — 

= 2i(xfcZ„+i - ZkXn+i). 

The last equality follows because the sum is a telescopic sum. D 

Since 1 = id is independent oft, PropositionfTil gives the derivative of the correction 
(^7r("))2_7^(n) |-Q g^ projector. As above, this gives an easy way for estimating the correction 
itself provided we have some clue how to choose the constant of integration. 

Our next task will be to construct, from the superadiabatic projections 7r'^"\ the 
transformation T" into the n**^ superadiabatic representation. By ([M)) and ((89| . tt*^") is 
almost a projection. Thus it has two eigenvectors w„ and Wn- Let 



be the eigenvectors of ttq, and write 

Vn = ckfo + /3wo, Wn — awQ — j3vQ (a, /3 G C). (109) 

We make this representation unique by requiring that < a G E and a^ + |/3p = 1, 
where the latter condition ensures that the determinant of T" is constant equal to 2, 
which is the value for Tq. 

Let Te be the linear operator taking (u„,w„) to the standard basis (61,62) of R^ , 
i.e. _ 

/ a+P a-0 \ 

T-i = (w„,u;„)= ^ Vp Vp (110) 

and _ 

where all vectors are column vectors. T^ diagonalizes ■k^"-\ thus 

T,7r(")T-i = ^ ^ Pr,' M . (112) 



A 



where Ai,2 are the eigenvalues of tt^"'. Although a,/3 and Ai^2 depend on n, e and f, we 
suppress this from the notation. 

Lemma 4. 

T,r-=(; f), ™. r,(T-.y = (;: f ) 4^' ( f .^ 
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Proof. Just compute, using 




To = t:^^ ^\ (113) 



and 



/ a'+l3' a'-l3 \ / / _ a+P a~0 

(T"^)' = Vp VP_ +£_ VP Vp 



D 



It will turn out that (3, a' a, and /3' are small quantities, A']^,A2, and A2 are even 
much smaller, while a^ and Ai are large, i.e. of order 1 + 0{e). This motivates the form 
in which we present the following result. 

Proposition 15. Suppose Ai 7^ A2. Then for each n e N, 

V A1-A2 (^n+l+Zn+l) -2 y 

^ ^ / |/?|2 - e(Im(^/3') + i0'aRe/3) i3(a:«+i + ^n+i) \ 

'"^'^(a;„+i-z„+i) -|/3|2+e(Im(^/3')+i0'aRe/3) y 



Ai-A: 



Proo/. Let us write T^iied^ - ff)^^"^ = (Mij), i,j G {1,2}. M14 and A/2,2 are calcu- 
lated in a straightforward manner, using Lemma |4] together with the fact TqHTq = 



diag(i,-i): 
















T.iied^ 


- H)Tr' 


= ied^ + ieT,T^'ToiT~^y - T.T^'ToHT^^ToT-' 








= lea^ +ie „ 

\ -P 


a J\ P'- 


-\e'p 




\&P 








I f a P^ 
2{-P a , 


1,0 -1 j 


if" 


f) 




Carrying out the matrix multiplication yields 














Mi,i = 


ied^ + ie{aa' + PP' 


- e'aReP) - 


\{a^- 


- i/3n 




Tl 


TIC rriTT-OQ fVid 


A/2,2 = 

rliarrriTial n 


ied^ + ie{aa' + p'^ 


+ e'aReP) + 


. Ifl|2 _ 


-|/3|^). 

1 cr, t 


Vicif rM 



2 



1 — 2|/3p and ad = —\{PP +PP'). Although we could get expressions for the off-diagonal 
coefficients by the same method, these would not be useful later on. Instead we use (J112p , 
i.e. T-'^D = 7r(")T-i together with ([TOTl) and obtain 

T,{ied^ - H)T-'D = DT.iied^ - H)T-^ - e"+ir,(z„+iX + x,,+^Z)T-\ (114) 
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By multiplying (|114p with e^e* from the left and by e^e^ from the right, j,k S {1,2}, 
J ^ 0, and rearranging, we obtain 

(Afc - A,) e, e* T^ed^ - H) Tf^ek el = (115) 

= -£"+^ej e* Te(z„+iX + x„+iZ) T^^Ck el . 

From the equalities TqXTq^ " ( - 1 V )' ^"^^o"^ " ( - 1 ) ^^'^ Lemma H we 
obtain 

' ^ I -(«^-/3^) «(/3+^) r 



T.ZT-^ = 



-1 



a{l3 -13) a2 + /?■ 



2 . ^2 



-(a2+/32) -a(/3-/3) ^ ' 

The expressions for Mi ^2 and A'h^i now follow from pisp . D 

We now use our results from the previous section to express a,/? and A1.2 in terms 
of Xk,yk and Zk, k < n. Let us define 

X ^ xKe,0 = ELi£'2;fe(0, (116) 

n EE ?7Ke,0 = ELi£'yfe(0, (ii7) 

Moreover, let 

g ^ g{n, e, = ELi£"+'5«+i,fe(C) (119) 

be the quantity appearing in (jlOSp . 

Lemma 5. The eigenvalues o/tt^"' solve the quadratic equation 

Ai_2 - Ai,2 -.9 = 0. 
Proof. By ()112p and Proposition [14] we obtain 

^'-^1 ^/^J=r.((7r("))^-.("))T-i=T,glT-i=(^ -^ ^ 

D 
Lemma 6. 

a^(Ai - A2) = 1 - ?/ - A2, and a/3(Ai - A2) = -X ~ C- 

Proof. From (jl09p and a^ + \f3\'^ = 1 one obtains 

iiQ = av„ - (3wn ■ (120) 
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Thus, using again (|109p . we find 

tt^'^^vq = Xiavn - X2Pwn ^ (Xia^ + X2\P\'^)vo + {\i - \2)a(3wo 
= (a^(Ai - A2) + A2)wo + (Ai - A2)a/3wo- 

On the other hand, from ([M)l and (ISS)) we have 



^("^=7ro + ^£'^(xfeX + yfey + ZfcZ), (121) 

fe=i 

and since Xvq = Zwg = — wq, ttoWo = wq and Yvq = —vq, we find 

7r(")wo = (1 - 77)^0 - (X + C)wo. 

Comparing coefficients finishes the proof. D 

We can now prove Theorem [U which we restate here in shghtly diflFerent form. To 
make the connection, note that a;„ is purely imaginary and Zn is purely real, so that 

indeed Xn + Zn = ~Xn + ^n- 



Theorem 7. Let en > Q be sufficiently small. For e G (0, eq] assume there is a bounded 
function q onM. such that x(^), r]{^), C(^) and their derivatives x' (0 ^ v' iO t C iO '^'"^ "^^^ 
bounded in norm by eq{^). Then 

TSed^ - H)T-^ = (122) 

/ k + 0{e^q) -£"+i(a;„+i-^„+i)(l + 0(e<z)) \ 

" '^ * ~ V £"+H^«+i + zn+i) (1 + 0{eq)) -\ + 0{e^q) ) ■ 

Proof. From p2ip and our assumptions it follows that tt*^"^ — ttq = 0{eq). Thus Ai = 
1 + 0{£q) and A2 — 0{eq), and from Lemma[n]we infer g = 0{eq) and 

Ai = i (1 + v/r+I^) , A2 = i (1 - v/TT4^) . 

Since Ai — A2 7^ 0, Lemma IH] yields 

2 l + ^/^+4^-2r; -^ - C 



2V1 + 4g VI + 45a 

Hence a^ = 1 + ©(eg), and /3,/9' and aa' — (a^)'/2 are all 0{eq). We now plug these 
into the matrix R in Proposition [151 This shows the claim. D 

6 Asymptotic behaviour of superadiabatic represen- 
tations 

The next step in the program is to understand the asymptotic behavior of the off-diagonal 
elements of the effective Hamiltonian in the n"^ superadiabatic basis for large n. Accord- 
ing to (|122|) this amounts to the asymptotics of Xn and z„ as given by the recursion 
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from Proposition [131 It is clear that the function 9' alone determines the behavior of this 
recursion. As in jBeTeli IBeTe2j it will in fact be the poles of this function closest to the 
real axis that play the dominant role; thus we first investigate the special case 

0'(O=7f— V + ^^) . (123) 

and then turn to a perturbation as allowed in Assumption 2. Fortunately, all the changes 
relative to [BeTeli IBeTe2) occur in the error terms of the explicit part stemming from 
(|123p . and so we have to redo only a small part of the work. 

We use (I106P in order to determine the asymptotics of z„. From Proposition [T3l 
together with ((97)) - ((99)) it is clear that y„ must go to zero as ^ ^ ±oo. This fixes the 
constant of integration in (|106p , and we arrive at the linear two-step recursion 

^„+2(0--^ U{0-e'{0fj{s)z^{s)ds\ . (124) 



Let us define 



?' = -:^(/-/) with /(O- '^' 



ic"' e + i^c 

and let us put ctj = 1 — 2(5oj-, i.e. aj = — 1 for j = and (Jj = 1 otherwise. We obtain 
Lemma 7. For each m > 1, 

m— 2 / -, \ fc+1 



rn — 1 

= -f E'^^-2--'Im(r+i-^) (126) 

Im(/™)' - _^Re(r+i), (127) 

Re(/")' = !!^lm(/'"+i). (128) 



Proof. We have f f — \{f + f), and thus 



and 



/V = ^/'-'(/ + /) = ^(/' + /'-V) 



/ ra — l /I \ k-\-l / 1 \ '"" N 



07" = ^ (/'"+' - /"7) = ^ ( /'"+' - E ( 2 j ^""' " ( 2 j ^ ) ■ (^^^) 

Taking the the real part of (fT29)) yields p26)) . To prove p27)) and p28)) . it suffices to 
observe that (/'=)' = i|:/''+^). D 
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"o 


- 1, af^-0, 




af+^^ - 


(n + l)n \^^ -^^ ^ 


i fc 
' ^^ n-A: ^ 

fc=0 m=0 



Proposition 16. For each even ti G N and j = 0, . . . ,n — 1, let the numbers a" be 
recursively defined through 

(130) 
afc-™a(;'M (131) 

^ ' \ fc=0 m=0 / 

for j < n and a^^'^^ — . Then 

2„ = -7^^^-^5]2-^-4")Re(r-^) {neven), (132) 

y„ = -7'^^^7^E2"M^E^^-^-«i"MRe(r"') (r^ e«en) , (133) 
Proof. We proceed by induction. We have xi = i9' /2 ~ P^Ini(/), and thus by (|100p and 

dnsD, 

Z2 = -ix'i = --^Re(/2). 

Sc 

This proves p32p for n = 2. Now suppose that (|132p holds for some even n € N. Then 
by (|100p and (|128p , (|134p holds for n — 1 . To prove (|133p for the given n, we want to use 
(jlOip . (|126p and the induction hypothesis on z„ yield 

0'.„ = __^2(!L_i):^^(n) g ^^2-('=+^-)lm(r+i-0+'=)) = 

= -7^^^E2~" E^™-.4"MMr+^-'"). (135) 

Since (|135p only contains second or higher order powers of /, it is easy to integrate using 
(fT28l) . Let us write 



6W^^^,_^.„(«). (136) 

3=0 

Then by p^ we obtain 

0'(s)z„(s)ds = -7'^^^^ E 2-"6L"^Re(r-"), 

-OO SC ™_n 
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proving (|133p for n. It remains to prove (|132p for n + 2. We want to use (|124p . and 
therefore we empfoy (I126P and our above calculations in order to get 

e'{t) f 0'(s)z„(s)ds- 

•/ — oc 

/' _ 1M ^~^ n-j + 1 



BydEzD, 



^^ j=0 \k=0 






Now we take the difference of the last two expressions, take the derivative of that and 
obtain 



Z„+2 = -7- 



(n-1)! 



n+2 



^ 2-^in +l-j)i in j)a^ - 7' ^ a.-^fc^ Re(r+2-^). 



j=o \ fc=0 / 

Comparing coefficients, this proves (I132p for n + 2. D 

We now investigate the behavior of the coefficients a, as n -^ 00. 

Proposition 17. Let a^" be defined as in Provosition UR 

(b) There exists Ci > such that for all n ^N 

\a^"^\<C, 



(n)| - „ Inn 



n-1 
(c) For each p > I there exists C2 > such that for all 7i G N 



supp ^|a,"^| < 



C2 
j>2 " n-1 

(2) 



Proo/. (a) By (fT30l) . a^ ' = 1, and 

(«+2) _ („) (^1 _ 2! 

Comparing with the product representation of the sine function ( [AbSt| . 4.3.89) 

00 / 
sin(7ra;) — ttx I I 1 ^ 
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we arrive at (a). 

(b) Put a„ = (n - l)a^"\ Then by p3T|) . 






a,i+2 = an 1 



thus for n — 1 > 7, we have 



|an+2| < |q;„| +7 - 



7^croCTi 



1 1 



1 1 



n n — 1 



,(») 



n n — 1 / mG 



max a, 



(m), 
I' 



which shows (b). 

(c) Put c^" — [n — l)p^%^" , and c^"^ = inaxj>2 \cj\. We wiU show that the sequence 

c*^"^ is bounded. We have 



^(n+2) 



11-/ -J' '= 

n(n - 1) I ^ ^^ ■' ^ 2^ n - fc ^ '""^ " 

^ ' \ fe=2 m=2 

\ fe=0 fc=i / / 



Now 



and 



fe=2 m=2 



<c(")- 



n~Jip~ 1)^ 



p-.^^_<ll±l>2< 



1 



fc=0 



n — k n — j (n — j)lnp 



(137) 



(138) 



(139) 



We plug these resuhs into (|137[) and obtain 

,(n+2)| ^ j„) / (" + l-i)("-j) I (n+l~j)7V 1 
'' ~ \ n{n — 1) (n — j)(p — 1)^ n{n — 1] 



1 (n + l-j>^7' nJn), ^ ,Jn) 
n [n — j){p — 1)^ Inp 



2,„ji«ri + i«ri)- 



By (a) and (b), Aq and a]^" are bounded. Taking the supremum over j > 2 above, we 
see that there exist constants Bi and B2 with 



^(„+2) < ^(n) 



n-2 



Bi 



n(ri — 1) 



B2 



hence 



,("+2) _ ,(n) < i ("("-2 + -^^ c(") + B, 



Now let n — 1 > Bi. Then for c^"-' > i?2, the above inequality shows c*^"+^-' < c'-"-', while 
for c(") < B2, c("+2) < c(") + B2/n < ^2(1 + 1/n). Thus c(") is a bounded sequence. D 
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Corollary 4. Let b " be given by Iil36\) . Then for each p > 1, there exists C3 > such 
that 

-jAn) ^ Gz 

supp ■'b, < -. 

j>0 ■' n-l 

Proof. For j < n — 1, we have n — 1 < j{n — j), and thus Proposition [17] (c) gives 
-ihi") ^ P^^ ( { (") I (")\ , ^2 P^ , ,_i -.A , 



n— 1 p— l?i — 1 71 — 1 



D 



Given the results of PropositionfTTI and Corollary |4] together with Lenima[7j we can 
now derive Theorems [2] and [3l The proofs are word for word the same as in [BeTelj and 
|BeTe2] from here on, and so we omit them. 
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